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Time-dependent electron transport through a quantum dot and 
double quantum dot systems in the presence of polychromatic 
external periodic quantum dot energy-level modulations is studied 
within the time evolution operator method for a tight-binding 
Hamiltonian. Analytical relations for the dc-current flowing 
through the system and the charge accumulated on a quantum 
dot are obtained for the zero-temperature limit. It is shown that 
in the presence of periodic perturbations the sideband peaks of 
the transmission are related to combination frequencies of the 
applied modulations. For a double quantum dot system under the 
influence of polychromatic perturbations the quantum pump effect 
is studied in the absence of source-drain and static bias voltages. 
In the presence of spatial symmetry the charge is pumped through 
the system due to broken generalized parity symmetry. 



1. Introduction 

Recently, considerable progress has been achieved 
in fabricating low dimensional systems and many 
experimental and theoretical works have been put 
forward. Especially interesting are quantum systems 
under the influence of external radio or microwave 
electromagnetic radiation perturbations where many 
interesting effects are observed like photon assisted 
tunneling (PAT) [1, 2], turnstile effects and photon- 
electron quantum pumps [3-5], conductance oscillations 
[6], and alike [7]. 

The symmetry of quantum dot (QD) systems (with 
no source-drain voltage) plays the crucial role as 
concerns electron pumping effects. Generally, one can 
consider symmetries like time-reversal symmetry, time- 
reversal parity and generalized parity [7]. 

A single electron pump based on asymmetrical 
couplings between QD and the left and right electrodes 
was considered in Ref. [4]. The couplings were switched 
on and off alternatively from zero to maximal values (by 



means of additional electrodes) and it led to electron 
pumping. A similar effect can be achieved for dipole 
driving forces applied to a double QD system (in the 
large gate voltage regime) or to quantum wires. In 
this case one QD site is driven by the external dipole 
interaction which is out of phase in comparison with 
the perturbation applied to the second QD site (the QD 
sites are not driven in homogeneous way), e.g. [7-12]. 
However, in the presence of spatial symmetry (and in 
absence of a source drain and static gate voltages) it 
is also possible to pump electrons but the symmetry 
must be broken in a dynamical way. The easiest way 
to break the time-reversal symmetry is to add a second 
harmonic to the driving system; i.e. a so called "harmonic 
mixing" drive [7,13,14], or in general the second external 
perturbation with an arbitrary frequency, [15]. In such 
a case, depending on the parameters of these two time- 
dependent perturbations perturbations the generalized 
parity can be broken and nonvanishing current can flow 
through the system [7, 13, 16]. 

There are few studies which address the electron 
transport through low-dimensional systems in the 
presence of several polychromatic external perturbations 
with arbitrary frequencies. Due to numerical problems, 
most of them concentrate on the case of commensurate 
frequencies or only bi-harmonic perturbations. It was 
shown that the external bi-harmonic time-dependent 
perturbations can be used to control the noise level 
in quantum systems [9, 17] or, as well, for routing 
optically induced currents [18, 19]. The shot noise for 
a single level quantum dot under the influence of 
two ac external perturbations (coherent or incoherent) 
was analyzed in Ref. [15]. The coherent destruction of 
tunneling [20] and the associated dynamical localization 
in quantum dots under the influence of a time-dependent 
perturbation with many harmonics were investigated 
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in Ref. [21]. Moreover, dissipative quantum transport 
in one or two dimensional periodic systems that are 
subjected to electric harmonic mixing perturbations (bi- 
harmonic) were studied in Refs. [22-24]. The nonlinear 
signal consisted of e.g. two rectangular-like driving forces 
which in turn allow to control overdamped transport in 
Brownian motor devices, [16,25-27]. 

In this paper we shall investigate the influence of 
polychromatic time-dependent energetic perturbations 
with arbitrary (commensurate and incommensurate) 
frequencies applied to the QD or double QD system 
attached to leads for charge accumulated on the QD 
and the time-averaged, dc-current flowing through 
the device. For a double QD system we propose a 
quantum pump which is based on a scheme which 
mimics closely a dipole-like perturbation. Thus, our work 
can be treated as a generalization of the studies of 
the electron transport through a QD or double QD 
systems affected by one external perturbation or bi- 
harmonic electric time-dependent ac-perturbations with 
arbitrary frequencies. A tight-binding Newns- Anderson 
Hamiltonian and evolution operator method are used in 
our calculations. 

The paper is organized as follows. In Sec.[2]the model 
Hamiltonian and theoretical description of a single level 
quantum dot are presented. Also analytical relations 
for the time-averaged current and time-averaged charge 
on the QD are obtained and numerical results are 
depicted and interpreted. In Sec. [3] the current through 
a double QD system is obtained and the pumping effect 
is discussed. The last section, Sec. [4] is devoted to 
conclusions. 

2. Single-Level Quantum Dot 

2.1 Theoretical description 

In this section, starting from the second quantization 
Hamiltonian and using the evolution operator method 
we obtain the charge accumulated on a QD and the 
current flowing through the system under the influence 
of many external time-dependent perturbations. Our 
system consists of a single level quantum dot and two 
connecting electron electrodes, left (L) and right (R). 
The total Hamiltonian is then given by H = Hq + V, 
where 

H o= £ k a c L c k a + £ d(t)c+c d , (1) 

ka=L,R 

V = H + E V %R C % R C « + h - c - ( 2 ) 

kL kR 



The operators c^ a (ci ), Cd(c^) are the annihilation 
(creation) operators of the electron in the lead a (a = 
L,R) and at QD, respectively. The QD is coupled 
symmetrically to the leads through the tunneling 
barriers with the transfer-matrix elements V^ L and 
Vt R (hopping integrals). For the role of asymmetric 
lead-"molecule" coupling see in Refs. [28, 29]. The 
electron-electron Coulomb interaction is neglected in our 
calculation, cf. [7,15,30]. 

External perturbations are applied to the QD (the 
QD energy level is driven in time by time-dependent ac- 
voltages). We consider here a harmonic modulation of 
the external energy level perturbations applied to the 
QD, i.e. 

n 

£d(t) =£d + ^A l cos(w i t + i ) (3) 

i=l 

where u>i , A j and </>j are the frequency, driving amplitude 
and phase of i-th perturbation. 

The current flowing through the system and charge 
localized at QD can be described in terms of the time 
evolution operator U(t, to) given by the equation of 
motion (in the interaction representation, h = 1), i.e., 

i^-U(t,t ) = V(t)U(t,t ), (4) 

where V(t) = U (t,t )V(t)U£ (t,t ) and U (t,t ) = 
Texp(ift o dt'H (t')y The knowledge of the 
appropriate matrix elements of the evolution operator 
U(t, to) allows us to find the charge accumulated on the 
QD, n d (t) (cf. [6,31,32]), which is given by: 

n d (t) =Y,Mto)\U d At,to)\ 2 , (5) 



where np(to) represents the initial filling of the 
corresponding single-particle states (/? = d, fcL, kR). The 
current flowing e.g. from the left lead can be obtained 
from the time derivative of the total number of electrons 
in the left lead, cf. [33]: 

j L (t) = -edn L (t)/dt, (6) 
where riL(t) can be expressed as follows: 

n L (t) =^n fc (t) =EE n "WI^(*' t »)| J ■ ( ? ) 

kL kL I 3 

Using Eq. [4] the following differential equations for 
Ud,p(t,to) and U% L g(t,to) matrix elements which are 



2 



ISSN 0503-1265. Ynp. $ia. atcypH. 2003. T. 48, N 9 



DYNAMICALLY BROKEN SYMMETRY 



needed to obtain the current and QD charge, can be 
written in the form: 



d 



^UdA t M= E Vd,a{t)Ua,d(f,to) 



ka=L,R 



^l^Mo) = W^(Mo) 

where nonzero elements of the function V are 

ft 



(8) 



(9) 



V it* = v L, d = V L ex p [* J to MO - *aXJ (io) 

Next, using the wide band limit approximation 
r Q (e) = 2^£fc Q V %o V?J{e - e ia ) = T a and assuming 
r one can find the following relation for 
Ud,d(t) matrix element (to = 0): Ud.d{t) = cxp (— Tt), 
and similar for U d % a (t) and U^ L ^(t): 



^kLJ3 



(*) 



iexp(-rt) / dt'V d j a (t')exp(rt') (11) 



(12) 



1 / dt ' V kL J?)Ud,f>(f) 



Note, that for t — ► oo the element Ud,d{t) tends to 
zero and thus the charge accumulated on the QD 
does not depend on the initial QD occupation np(to) 
for large t. Also the current flowing through the 
system is independent on nd{to) because the element 

3t2fc£ n «i(*o)l^fc£,d(*)| 2 tends to zero for 1 ~* 00 • 

Finally, the QD charge can be written in the form: 



n d (t) 



E n kM\u drka (t)\ 



(13) 



ka=L,R 



and the current through the system reads (e = 1): 



the QD-energy level, Eq. [3] and assuming 4>i — 0, the 
elements V d s , Eq. [10] can be expressed as follows: 



i—l rrii 



(15) 



and the solution of Eq. [TTJ for the evolution operator 
elements can be written in the form: 



-V* e** 6 * -6 **'* x 

ka 



E-E 



mi m-n 



J, 



»i (wi ) — Jrrin ( w") ' 

£rf — e? + — iT 

u ka 



(16) 



where mi G (— oo, oo) and is an integer number, fl = 
m\ui\ + ... + m n uj n and Ji is the Bessel function. 
Next, we obtain the dc current through the system, 
jo = (j(t)) = limr^oo *f; j-T/2 which can be 

symmetrized, j = O'l(O) = - A and 

finally we find: 



io - ^ f° (/*(e)-/z(e))T(e) 



(17) 



where Jl/r^) is the Fermi function of L/i? electrode 
and the transmission reads 



(18) 



T (e) = 1-E-EE-E'V, .•' - 

^'"1 (wi") •■■^ m « (uf) ^mj (t3t) '"^"^ (l3r) 

(e d - S + ft) 2 + r 2 

where 51 = m x u)\ + ... + m n u> n . Note, that 
for incommensurate frequencies the Kronecker Delta 
function S n _ n ' is nonzero only for m i — mi and thus 
the transmission can be written in the following short 
form: 



JL 



(t) = rnd(t)+Im\^n u (t )Vg L /t)U dM (t)) () = 2 v v J ^ (%) - J L (%) 
\& / 1 ' (e d -e + «) 2 + r 2 



(19) 



(14) 

Eq.[T3]and Eq.[l4]are very general relations which should 
be analyzed using Eq. \W\ and Eq. [TTJ The relation for 
the current, Eq. [14l has the structure which can be 
written by means of the transmission Tlr and Trl, 
i-e. 3L(t) = Y,U n %L$)TLR{t) - Y,% R n %R {Q)T RL {t), 
but note that in general TL R (t) ^ Trl(<). In order to 
obtain the current one should know the exact form of 
V d £ q function. Using the time-dependence relation for 



This equation is valid also for the case when there is 
large difference between frequencies u>i. Note, that the 
transmission, (Eqs. [18] and [19]), corroborates with the 
results obtained by means of Green's function method 
for a quantum wire driven by homogeneous external 
perturbations, [30]. 

The relation for the current, Eq.[l7l has the structure 
of the Landauer formula. We note that the structure 
in Eq. (17) involves, in clear contrast to a electric 
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field dipole perturbation [7] , no inelastic photon-assisted 
tunneling events. This is so because with the time- 
dependent energy level perturbation used here the long 
time-average of Tl_r(£) and TRL{t) equals TLn(e) = 
Trl(c)- For zero temperature and for incommensurate 
frequencies of external perturbations the dc-current can 
be obtained analytically from the relation: 

mi m„ v ' \ n / 

( £ d -/ij,-fi £ d -fi R -ft\ 
arctan arctan 

V r v ) 

Note that for = M-R the dc current is zero. 

The corresponding analytical relation for the 
accumulated quantum dot dc-charge; i.e. no = (nd(t)) = 

limT^oo <f J T j 2 /2 n d(t')dt' reads: 

- - s (£)■•-*•(£)'< < 21 > 

mi m„ \ 1 / V n / 

( £d - Ml - n e d - Mr - n \ 
I 7r — arctan — arctan — I 

It is worth noting that for Ed << fi (ed » (i) the 
charge accumulated on the QD is maximal (minimal). 
The relations for the dc current, Eq.l20| and for the QD 
charge, Eq. [2TI constitute the main analytical relations 
of this section. 

2.2 QD accumulated charge and dc current 

In this section we analyze the QD charge and the 
dc current flowing through a quantum dot driven by 
polychromatic perturbations. All energies are expressed 
in the units of T° and in order to obtain rather narrow 
sidebands peaks we assume T = 0.2T (taking the 
unit of energy T° — 0.05eV it corresponds to T = 
O.OleF). For larger T all sideband peaks are wider and 
it is difficult to observe many-perturbation effects. The 
current and the conductance are given in units of 2eT°/h 
and 2e 2 /h, respectively. Moreover, we show numerical 
calculations for two external perturbations case but the 
generalization for more perturbations is obvious. 

In Fig. 1 the QD charge (upper panel) and the dc 
current flowing through the system (lower panel) are 
shown for two external perturbations applied to the 
system (cji = 3, CJ2 = 8, Ai = 4, A2 = 8) - thick lines. 
The frequencies are commensurate but there is large 
difference between them and the transmission obtained 
form Eq. [18] and Eq. [19] are almost the same. Physically 



it means that eight-photon adsorption/emission process 
based on the third sideband peak should occur to play 
the role in the transmission, (8ui\ = 3cj2), but this is 
unlikely process. The broken lines correspond to the 
single external perturbation case, i.e. A2 = 8 (Ai = 0) 
- thin broken lines and Ai = 4 (A2 = 0) - thick broken 
lines, respectively. The chemical potentials are rather 
small, fiL = — (iR = 0.1, and the current peak for ed = 
(lower panel) is observed for mono and polychromatic 
cases (it appears also for the time-independent case). 



0.8 




-0.2 I-" ' 1 ■ "-i— ■ — l 



-2 2 4 6 8 10 12 

e d 

Fig. 1. QD charge (upper panel) and the dc current (lower panel) 
as a function of ed for Ai = 4 and A2 = 8 (thick lines), Ai = 4, 
A2 = (thick broken lines) and Ai = and A2 = 8 (thin broken 
lines), respectively. The other parameters are &i = 3, u>2 = 8, 
fiL = —fJ.R = 0.1, r = 0.2. The thick (thin) broken line on the 
lower panel is shifted by -0.1 (-0.2) for better visualization. 

For only one external perturbation applied to the QD, 
lo = 3 (or w = 8), the sidebands peaks are visible for 
ed = iku>, where k is an integer number. However, for 
the case when two external perturbations are applied 
simultaneously to the QD, additional sideband peaks 
appear. Of course the single peaks from the first and 
the second fields are still visible i.e. for Ed = uii, 2u>i or 
u>2- In Fig. 1, additional dc current peaks are indicated 
by points A\, A2 (first-order sidebands) and B\ (second- 
order sideband). A\ (A2,B\) sideband peak appears for 
£d = ^2 - uji (ed — Lo 2 + wi, Ed = UJ2 — 2cji) and 
is related to the peak for ed — 1^2 — 8 (but not to 
the main peak for ed = 0). It is worth noting that 
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although the frequencies u>i and U2 are not equal to 5 
or 11 we observe the sideband peaks for these values 
of Ed- In general we observe sideband peaks for Ed — 
±kiUJi ± /C2W2, where k\ and ki are integer numbers. 
The structure of the dc current curves are reflected also 
in the charge accumulated on the QD (upper panel). 
The charge decreases with Ed but there are many steps 
which are related to the current sideband peaks. Also 
additional sidebands i.e. points A\, A<x and B\ from the 
lower panel, are visible on the charge curve (Fig. 1, upper 
panel, thick line). 



0.04 



0.02 




4 6 
A-|, A 2 =10-A 1 

Fig. 2. The dc current as a function of the amplitude Ai (A2 = 
10 — Ai) for ed = 5 (thick line) and Ed = 10 (thin line), ui\ = 3, 
u>2 = 8. The thin (thick) broken line corresponds to one external 
perturbation case Ai, lu = 3 with A2 = (A2, lu = 8 with Ai = 0) 
and Ed = 5; hl = —)J-R = 0.1, T = 0.2. 



conclude that one external perturbation slightly changes 
the current in this case. It confirms that the current 
maximum, which appears for Sd = 5 (thick solid curve), 
is due to a combination of two external perturbations 
applied to the system. The similar conclusions are valid 
for other £d = ±fciWi ± fc2^2 (&i,2 7^ 0). 



0.15 



0.1 



0.05 



i.15 
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Next, in Fig. 2 we study the dc current as a function 
of the driving strengths (amplitudes) of two external 
perturbations applied to the system. For the second 
signal the driving strength decreases with the amplitude 
of the first perturbation, i.e. A2 = 10 — Ai. The thick 
(thin) solid line corresponds to ej, = 5 (54 = 10). As 
one can see for Ed = 10 the dc current is very small 
and almost independent on the amplitudes Ai and A2. 
This conclusion is valid for every Ed which is not any 
combination of u)\ and cl>2. For Ed = 5 (Ed = u>2 — irfi) 
the current is minimal for Ai = and Ai = 10 but 
for the mixed regime of Ai and A2 it is characterized 
by a local maximum, cf. thick solid line. This is a very 
interesting effect because one can control the current 
flowing through the system by applying additional time- 
dependent perturbation. Note, that the maximal value 
of the current is a few times larger than for the case of 
one external perturbation, i.e. for Ai = or Ai = 10. 
For only one external perturbation applied to the system 
the dc current is shown by the broken lines - the thin line 
for A2 = as a function of Ai (u> = 3) and the thick one 
for Ai = as a function of A2 (u> = 8), 64 = 5. One can 



Fig. 3. (Color online) The dc current as a function of Ed and LU2 
for Ai = 0, A2 = 8 (one external perturbation - upper panel) and 
Ai = 3, A2 = 8 (two external perturbations - lower panel). The 
other parameters are: oji = 3, = — Hr = 0.1, T = 0.2. 

To analyze the structure of sidebands peaks for 
polychromatic perturbations in Fig. 3 we show the dc 
current as a function of the QD energy level, Ed, and the 
frequency of external perturbation L02 for Ai = (only 
one perturbation applied to the system with A2 = 8) 
- upper panel, and for two external perturbations with 
Ai = 3, u)\ = 3, A 2 = 8 - lower panel. For one, as well 
as for two time-dependent perturbations applied to the 
system we can distinguish between the adiabatic regime 
for L02 < 1 and the non-adiabatic one for W2 > I, For 
L02 > 1 one can observe the main current peak for Ed = 
and sidebands which are visible for Ed = ±ku>2- The 
situation is more complex for two external perturbations 
applied simultaneously to the QD (Ai = 3, u)\ = 3, 
A2 = 8). One observes that each light line from the 
upper panel possesses satellite lines which are localized 
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in the distance ±ku>i from the original lines. Thus we 
observe very reach structure of the dc current in the 
presence of polychromatic perturbations. 



2.3 Results: Transmission asymmetry and 
phase dependence 



0.15 



0.1 




0.05 - 



1 1 r 



to 2 /(OH =2 



C0 2 /CL)jk = 1 



to 2 /(Oi=0.5 





M | V ^ i ^ - ^ i 
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For one external perturbation applied to a QD as 
well as for two external perturbations with the same 
frequencies one observes fully symmetric current curves 
around the value Sd = 0. However, for two external 
time-dependent perturbations the dc current (or the 
transmission) can be asymmetrical which appears for 
the case of commensurate frequencies. In Fig. 4 we 
show the dc current obtained for the frequency ratio 
ui^/oji = 2, 1 and 0.5 using the transmission relations, 
Eq. HH1 (solid lines) and for the case where there is an 
infinitesimal shift of the frequency ratio from an integer 
number, Eq. [19] (broken lines), e.g. for UJ2/0J1 — 2 we set 
u/i = 3 and LO2 = 5.999. Note, that all broken curves 
in Fig. 4 are symmetrical. The solid lines, however, are 
asymmetrical around Sd = except for the case loi = u>2 . 
Notably, for oj% — uj 2 the system is equivalent to a single, 
effective external perturbation case with the effective 
driving amplitude, A = Ai + A2, and thus the current 
curves render themselves symmetrical again. The current 
obtained for incommensurate frequencies differs from the 
ones of commensurate frequencies. It is worth noting 
that the structure of the dc current curves for fixed u)i 
depends on the amplitudes of external perturbations but 
the positions of peaks remain unchanged. 



Fig. 4. The dc current as a function of £<j obtained for 
commensurate frequency ratios 1^2/(^1 =2,1 and 0.5 according 
to Eq. [17] with the transmission given by Eq. \TE\ (solid lines) and 
for a slightly off-commensurate frequency; i.e., uj\ = 3,^2 = 5.999, 
see Eq. [19] (broken lines). All parameters are the same as in Fig. 3. 

Next, it is necessary to explain the transmission 
and current asymmetries for the case of commensurate 
frequencies. In this case to obtain the transmission one 
should use Eq. [18] which is expressed by means of four 
Bessel functions, an energy dependent factor and the 
Kronecker Delta function. For commensurate frequencies 
this Delta function produces off-diagonal elements of 
the Bessel functions which are multiplied by the energy 
factor. The energy factor depends on the order of the 
Bessel functions and introduces asymmetry for positive 
and negative values of £d (it means that the energy 
factor possesses different values for positive and negative 
integer order of the Bessel functions). The off-diagonal 
elements of the Bessel functions are also responsible for 
different values of the dc current obtained for o^/^i = 1, 
according to Eq. [18] and Eq. [T9] It is worth noting that 
in our calculation we assume that the QD energy level is 
driven according to the formula with the cosine function 
(i.e. the cosine function is an odd function) and the 
external time-dependent perturbation satisfies the time- 
reversal symmetry, a(t) — a{—i). 
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Fig. 5. The transmission versus the energy for two external 
perturbations applied to the QD (cji = 3, lu2 = 6, Ai = 3, 
A2 = 8) with phases: rpi = <p2 = 0, ir/2, it (upper panel) and 
4>i = 0, it/A, it/2 and <p2 = (lower panel). 

In Fig. 5 we show the transmission for nonzero phase 
factors of the external perturbations i.e. 4>i and 2 , cf. 
Eq. 03 The upper panel shows the results for the case 
when both phases are changed in the same way i.e. 
4>i = 4>2- For 4>i — (thin solid line) the asymmetry 
in the transmission is observed, cf. also Fig. 4. If both 
phases are equal to ir/2 (broken line) the transmission 
is symmetrical (as a function of the energy) - in that 
case the external signals are sinuses. For <pi = tt 
the transmission is asymmetrical again. The period of 
the external perturbation for the case of the same 
phases, <j>-y — $2, is equal to 2n. The situation is 
somewhat different when the phase of one external 
perturbation is constant (4> 2 = - lower panel). Here, 
only for </>i = ir/4 the transmission is symmetrical but 
for 0i = and n/2 this function is asymmetrical. 
Note, that the above conclusions are valid for the 
case of commensurate frequencies. For incommensurate 
frequencies the transmission curves are symmetrical, cf. 
Fig. 4. Moreover, as one can see, the transmission curves 
in the upper and lower panels are the same, only the 
phases of external perturbations are different. It means 
that the same effect can be obtained by changing only 
one phase parameter instead of driving both phases 
simultaneously. 



3. Double quantum dot system and pumping 
effect 

The asymmetry effect in the transmission obtained 
in the previous section can be used to construct a 
single electron pump [12] based on a two-level fully 
symmetrical system with no source-drain and static bias 
voltages applied to the system. In this section we analyze 
the electron transport through a double quantum dot in 
the presence of external perturbations. The Hamiltonian 
of our system can be written as follows, in close analogy 
to Eq.CDand Eq.EJ i.e., 

#o= ]T efc - ct 4 Q C^+ ei (i) C +c 1 + e2 (t) C + C2 , (22) 

ka=L,R 

V = J2 ^4l ci + E v kn c t R C2 + v ^ c * + h - c - ( 23 ) 

kL kR 

where V12 is the tunnel coupling (hopping term) between 
two QD sites. As before, we assume that there are only 
two external perturbations applied to the system and one 
can write the following time dependence of the quantum 
dot energy levels: 

£1 (t)=e 1 +A l cos(wii) + A 2 cos(w 2 i) (24) 

e 2 (t) = e 2 + Ai cos(wi£ + 4>) + A 2 cos(w 2 i + 4>) (25) 

where 4> is the phase difference between the external 
perturbations applied to the first and the second QD 
sites which play a similar role as dipole forces in one 
external harmonic field case, cf. [7, 9, 10]. We stress, 
however, that here the driving is chosen uniform in the 
sense that the amplitudes strengths for the driving of 
the two dots are identical. Note, that the role of phase 
difference between the first and the second external 
perturbation was considered in [13] for bi-harmonic 
perturbations and here we concentrate on the phase 
difference between both QD sites. The main reason 
we omit the phase difference between both external 
perturbations is that in the presence of spatial symmetry 
(si = e 2 , pll = Ma; i-e., in absence of a static gate 
voltage and a source-drain voltage), it is impossible to 
pump electrons through the system. By introducing the 
phase difference according to Eq. [24] and Eq. [25] we 
change the symmetry at the first and the second QD 
sites and thus electrons can be pumped in the presence 
of external perturbations and spatial symmetry. The 
general formula for the time-dependent current flowing 
from the left electrode can be written as follows: 
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hit) 



E 

kL 



A- TV 



where Vjg L 1 (t) is denned by Eq. [15] and the evolution 
operator matrix elements satisfy the following set of 
differential equations (and similar for U 1( ^ % R ): 



difference effect as it can lead to the electron pumping 
An the system with no source-drain and static bias 
^L,i(^i,fci(*) Jvoltages. We take into consideration phase difference, 
(j), between QD sites and for <f) = it/2 it leads to a 
^2gjiipole-like parameterization (the oscillations are out of 
phase) . The similar two-level system under the influence 
of one external perturbation for the case of high bias 
voltage, e\ « £2, was investigated in Ref. [34] and in 
the presence of phase difference between both external 
perturbations and dipole driving forces in Ref. [13]. 



d_ 

at 



-iV %L e^-^)^h _ V -U irkL (t) (27) 



d_ 

at 



-iVvie^-^jUt-Mu^t) 



oU 2rkL (t) 



(28) 



where /1 = ^ sin(wii) + sin(w 2 i), f 2 = ^ sin(wi* + 
(f>) + sin(cJ2i + (j)). Note, that for (j> = 2nk, and the 
same on-site energies £\ = £2 = £<2, one can obtain the 
dc current and the transmission analytically. In this case 
the current satisfies the Landauer formula and the time- 
averaged left-right and right-left transmissions are equal. 
This result is due to the above mentioned chosen uniform 
driving strengths. For incommensurate frequencies of the 
external perturbations the formula for the transmission 
simplifies and can be written as follows: 



T(e) 



(29) 



J 2 



((e d - e + ft) 2 - f - V£) +T 2 {e d -e + ny 

where ft = uj\ra\ + W2W2- For the case of commensurate 
frequencies the formula for the transmission assumes no 
simple transparent from, but is similar to Eq. [TH Note, 
that also for (f> — and E\ ^ £2 it is possible to analyze 
the set of differential equations analytically, Eq. [27] 
and Eq. (28] however the solution is rather complex. 
Generally, i.e. for ^ and e\ 7^ £2 we solve the 
set of differential equations for £^(2) %l(K\ numerically 
then put the solution into the relation for the current, 
Eq. [26] and time-average it over the common period 
of external perturbations. Thus, this procedure can be 
applied only for the case of commensurate frequencies. 
In our calculations we concentrate mainly on the phase 



If only one external perturbation is applied to the 
double dot system, the phenomenon electron pumping 
does not occur for the case of spatial symmetry e\ = £2, 
[7]. For nonzero gate voltage and <j> = tt/2 (i.e. in 
the presence of the phase and spatial asymmetry) the 
current can flow through the system, Ref. [34], but if 
there is no phase difference between both QD sites, 
= 0, the current vanishes again. In order to analyze the 
role of polychromatic perturbations on the symmetrical 
system i.e. e\ = £2 = (no applied static gate voltages) 
we depict in Fig. 6 the dc current as a function of the 
driving amplitude, A2, for two external signals and for 
the phase difference between QD sites 4> = = 2ir (thick 
broken line), <f) — tt (thin broken line) and (j> — tt/2 (solid 
line). For A 2 = the current is zero, independent on 
the phase cj>; in this case only one external perturbation 
is acting and no symmetry breaking taking place. For 
two external perturbations, and for no phase difference 
between the first and the second QD sites, the current 
also does not flow, being equal to zero for all A 2 (thick 
broken line). However, for (j) £ (0,2n) the current flows 
although there is no voltages applied to the system. 
In that case, we realize a sort of a quantum pump 
which works only in the presence of two time-dependent 
perturbations, implying that generalized parity is broken 
in a dynamical way [7, 13]. Note, that depending on the 
phase difference between both QD sites the dc current 
can be positive or negative, cf. the broken and solid lines. 
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DYNAMICALLY BROKEN SYMMETRY 
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Fig. 6. The dc current flowing through the double QD symmetrical 
system under the influence of two external perturbations versus the 
amplitude A2 for the phase difference <f> = = 2w (thick broken 
line), <j> = 7r (thin broken line) and <j> = n/2 (solid line). The 
other parameters are u>i = 3, Ai = 3, ui2 = 6, V\2 = 1, T = 0.2, 
£1 = £2 = 0, m = fiR = 0. 

A remaining question is how the pump current 
depends on the phase difference between both QD sites. 
In Fig. 7 we depict the dc current flowing through a 
double QD symmetrical system (ei = £2 = 0) as a 
function of the (relative) phase 4> for three values of the 
driving strength of the second perturbation, i.e. A2 = 3 
(solid line), A2 = 6 (thick solid line) and A2 = 9 
(broken line). In case of zero phase difference, = 0, 2ir, 
the current is zero although two external perturbations 
are applied to our system, cf. also Fig. 6. For values 
of 4> different from 4> — 0, 27r a finite pump current 
flows through the double dot system: Depending on the 
relative phase it can be either positive or negative. 




-0.02 1 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 

*[«] 



Fig. 7. The dc current as a function of the phase difference (p 
for two external perturbations u\ = 3, Ai = 3, u>2 = 6 and for 
A2 = 3 (solid line), 6 (thick solid line) and 9 (broken line). The 
other parameters are: V12 =1,T = 0.2, £1 = £2 = 0, fi^ = = 0. 



It is also possible to stop the pumping current all 
together for specific values of (f>. Note, that for very large 
amplitudes of external perturbation, A2, the current 
takes on positive values for almost for all 0, cf. the broken 
line. 

4. Conclusions 

Time-dependent electron transport through a quantum 
dot and double quantum dot system in the presence of 
polychromatic perturbations has been studied within the 
evolution operator method. The QD has been coupled 
with two electrodes and external, time-dependent energy 
perturbations have been applied to the central region. 

Analytical relations for the dc current flowing 
through the system, Eq. [20] and the charge accumulated 
on a QD, Eq.[2lJ have been obtained for incommensurate 
external perturbations. Also the analytical relations 
for the transmission have been derived, Eq. [18] for 
commensurate frequencies and for incommensurate 
frequencies, Eq. [T9] It has been found that sideband 
peaks appear for Ed — Y^i=i ±fet<*^ where fcj is an integer 
number and n stands for the total number of external 
perturbations. For the case of two external perturbations 
this condition can be written as Ed = ikiu>x ± fc2<x>2 and 
e.g. for u>i = 3 and 0J2 — 8 additional peaks appear 
for Ed — 2,5,..., cf. Fig. 1. In the presence of external 
perturbations one can control the dc current by changing 
the amplitude strengths of the acting perturbations, cf. 
Fig. 2. In the presence of multiple external perturbations 
the current is characterized by satellite peaks, cf. Fig. 3 
for two acting energy perturbations. 

Moreover, the dc current obtained for commensurate 
frequencies, e.g. for bi-harmonic perturbations, strongly 
differs from that one obtained for slightly different 
from commensurate frequencies, cf. Fig. 4. Also the 
asymmetry in the transmission and the dc current versus 
the quantum dot energy level for commensurate 
frequencies case is detected, see in Fig. 5. This 
asymmetry appears only for commensurate external 
perturbations applied to the system and is related to a 
phase difference between time-dependent perturbations. 
For a double QD system the analytical formula for the 
transmission has been obtained for incommensurate 
frequencies, Eq. [29] Also a electron quantum pump 
based on a fully symmetrical double QD system has 
been proposed in absence of source-drain voltages and 
static bias voltages for which the pump current varies as 
a function of relative phase shift <f>, cf. Fig. 6 and Fig. 7. 
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